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Constraints on Dark Energy from Galaxy Cluster Gas Mass 
Fraction versus Redshift data 

Lado Samushia 1,2 ' 3 and Bharat Ratra 1 ' 4 
ABSTRACT 



We use the lAllen et al.l ( 120081 ) galaxy cluster gas mass fraction versus redshift 
data to constrain parameters of three different dark energy models: a cosmo- 
logical constant dominated one (ACDM); the XCDM parameterization of dark 
energy; and a slowly-rolling scalar field model with inverse-power-law potential 
energy density. (Instead of using the Monte Carlo Markov Chain method, when 
integrating over nuisance parameters we use an alternative method of introducing 
an auxiliary random variable.) The resulting constraints are consistent with, and 
typically more constraining than, those derived from other cosmological data. A 
time-independent cosmological constant is a good fit to the galaxy cluster data, 
but slowly evolving dark energy cannot yet be ruled out. 



Subject headings: cosmology: cosmological parameters - cosmology: observations 
- X-rays: clusters 



Introduction 



Observations over the last decade have established that the cosmological expansion is 
accelerating. In the context of general relativity this requires that the cosmological en- 



ergy budget be dominated by dark energy (for reviews see, e.g.. ICope 



and et al. 2006 



Padmanabhad 120061 : UshaJd 120071 : IUzanll2007l : iLinderl 120071 : iRatra fc Vogelevl 120081 ) . This hy- 
pothetical construct — dark energy — is an enigma. It is not yet clearly established whether 
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dark energy is Einstein's cosmological consta nt A (e.g., |Peebleall984), or whether it evolve s 



slowly in time and varies weakly in space (e.g. JPeebles fc Ratra 



1988 



Ratra fc Peebleslll988f ). 



While Type la supernova (SNIa) apparent magn itude measurements as a function of red - 



shift indicate accelerated cosmological expansion faiess et al.lll998l : IPerlmutter et al.lll999h. 



SNIa data are as yet unable to unambiguously constrain dark energy (see, e .g., 



200 71: iNesseris & Perivolaropoulosl 120071 : IShafielool 120071 : IZhang et all [2007; 



Ishida et al.ll2008l ). Future SNIa data will improve the constraints (e.g. JPodariu et al. 



Alam et al. 



Wu et a 



2008 



2001ah 



and could resolve some of the current differences in results from different SNIa data sets. 

The results of the SNIa test are confirmed by a test based on cosmic microwave back- 
ground (CMB) anisotropy data that must assume the cold dark matter (CDM) model for 
structure formation (see Peebles & Ratra 2003 and references therein for a discussion of 



verse having negligible spatial curvature (see, e.g., 


Podariu et al. 2001b: 


Durrer et al. 2003: 


Mukheriee et al. 


2003; 


Paee et al. 2003: SDergel et al. 2007: Doran et al 


[2007a), under the 


assumption that dark energy does not evolve in time (e.g., Wright 2006: 


Tegmark et al. 


2006: Zhao et al. 


2007 


Ichikawa & Takahashi 2007 


Wang & Mukherjee 


2007 


). In combina- 



tion with low measured non-relativistic matter density (Chen & Ratra 2003b and references 
therein), negligible spatial hypersurface curvature indicates the presence of dark energy. 

There are many different models of dark energy^ In this paper we consider three: stan- 
dard ACDM, the XCDM parameterization, and a slowly rolling scalar field dominated one 
(0CDM). In the ACDM model the late-time Uni verse is domin ated by a cosmological con- 



stant A with time-independent energy density p\ (IPeeblesI Il984j ) . In 0CDM the dark energy 



is a slowly rolling scalar field 0; in the model we consider the scalar field potent ial energy 



density V(4>) oc (f)~ a , where a is a nonnegative parameter ([Peebles fc Ratralll988l ). We also 
consider the XCDM parameterization of the dark energy equation of state; here dark energy 
is modeled as a fluid with an equation of state which relates the fluid pressure p = u x p to its 
energy density p where u x is a time-independen t negative p arameter. This approximation is 



inaccurate in the scalar field dominated epoch (jRatralll99ll ). In all three models, other con- 



tributors to the Universe's current energy budget include CDM and baryonic matter. In the 



: In this paper we assume general relativity is an adequate description of gravitation on cos mological scales 



For d is cussions o f accel e rated cosmologic a l exp a nsion in modified g r avity mode l s see, e.g., [M ovahcd et al 



2007), Fav et a 



Demianski et al 



(20071). lAmendola et all (|2007h . I Wang et all (|2007l ). IWan et all (|2007h . lWang! (|2007ah . and 
( 20071 ). We also assume that dark energy and dark m atter are uncoupled. For discussions 



of cou p led or unif i ed dar k energy and dark m atter models see, e.g ., Bonometto et al. ( 2006h . Wu fc Yu 



(2007), 



Guo et al 



(2007). Wei fc Zhand (120071). a nd 



Olivares et al 



(200. 



For other m o dels see, e.g. 



Barenboim fc Lvkkenl (I2006I). Ir3rax fc Martini d2007h. butta fc SorbcT(|2007l) . IWu et all |2007l ). brande et al 



(|2007h . and lNeupane fc Schereri (|2007h . 
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0CDM and XCDM cases spatial hypersurfaces are taken to be flat, while spatial curvature 
is treated as a cosmological parameter in the ACDM model we consider. 

It is important to confirm and strengthen the SNIa and CMB test results by using 
additional techniques. This will allow for consistency checks as well as possibly identi- 
fying systematic effects in a particular data set. Other promising curren t tests include 



the angular size of radio sources and quasars as a function of redshif t (e.g.. IChen fc Ratra 



2003al ; |Podariu et al 
lensing 



2007 



2003 



Chae et al 



Qguri et al. 



Da 



2004 



fc Diorgovski 



Alcaniz et al 



2005 



2005 



Dalv et al.lhoOTf). strong gravitational 



Fedeli fc Bartelmannl I2007L Lee fc Ng 



2008h . measurements of the Hubble parameter as a funct i on of redshift 



e.g., ISamushia fc Ratrakood : ISen fc Scherrer1l2008l : lLazkoz fc Maierottdl20071 ; IvVei fc Zhang 



2007; Samushia et al 



ments (e.g., 



Doran et al. 



2007 |), and large-scale structure baryon acoustic oscillation measure - 



2007bl:IParkinson et alil2007i:IPercival et alil2~00~7l: 



Lima et a 



For reviews of the observational situation see iKurek fc Szvdlowskvi J2008h . IWangj d2007bh . 



200 



3). 



and lLazkoz et al.l (120071 ). Current data favor dark e nergy that does not evolve, but do not 



yet strongly rule out evolving dark energy (see, e.g.. iRapetti et al.ll2005l ; IWilson et al.ll2006 
Davis et alil2007h . 



In this paper we use new galaxy cluster gas mass fraction versus redshift data (jAUen et al 



20081 . hereafter A08) to constrain parameters of t he three dark ene rgy models mentioned 
above. Earlier galaxy cluster gas mas s fraction data ( Allen et all2004) have be e n used to con- 
strain these and other models (e.g., Chen fc Ratral 120 04; Ra petti et al.ll2005l ; I Alcaniz et al. 



20051 : ICannata fc Kamenshchikibood : Izhanlhood : IWei fc Zhang! l2007h 



If rich galaxy clusters are large enough to have matter content that fairly samples that 
of the Universe then the baryonic to total mass ratio in clusters is equal to the cosmological 
baryonic mass fraction ratio of fib and fi m , the baryonic and nonrelativistic mass density 
parameters. The baryonic mass in clusters is dominated by the x-ray emitting gas which 
can be measured through x-ray observations. Combined with an estimate of fib from other 
observations, these measurements can be used to constrain fi m . The gas mass fraction in 
large relaxed clusters is expected to be independent of redshift and, since the observed x-ray 
temperature depends on the assumed distance to the cluster, this can be used to constrain 



Ettori et al.ll2006l : iNagai et al. 



cosmological parameters (ISasakil Il996l: |Penlll997 ; for r ecent discussions of the test see, e.g. 



20071 : ICrain et al. 



20071). 



In addition to depending on cosmological parameters, in a given model the predicted 
gas mass fraction depends on a number of "nuisance" parameters that have to be marginal- 
ized over to derive the probability distribution function for the cosmological parameters of 
interest. We note that the likelihood function depends on certain functions of the "nuisance" 
parameters and so introduce auxiliary random variables to describe these functions. This 
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technique helps us to significantly reduce the computational time. 

In next section we outline our computations. In Sec. 3 we present and discuss our results 
and conclude. 









_d A (z)_ 



2. Computation 

In our computations we follow A08. For a given cosmo logical model we compute pre- 
dicted values of the gas mass fraction 

s l + s (l + a s z) 

as a function of cluster redshift z, 4 cosmological parameters (h, Q m , fl^, and a parameter 
p, described below, that represents the dark energy model), and 7 parameters (sq, bo, a 
a&, K, A, 7) represented by Q and related to modeling the cluster gas mass fraction 
Here h is the Hubble parameter in units of 100 kms _1 Mpc _1 , d^(z) is the angular diameter 
distance computed in the reference spatially-flat ACDM model (with fl m = 0.3, ^a = 0.7, 
and h = 0.7, where Q\ is the cosmological constant density parameter), and dx(z) is the 
angular diameter distance computed in the cosmological model of interest and depends on 
the model and h. Since the angular diameter distance dx(z), and so f^,, depends on the 
assumed dark energy model, we can compare predicted values of the gas mass fraction with 
measurements for clusters at redshift Zi by constructing af = Si(/gas( z «) ~ fga$( z i)) 2 / (T i 
function (erj are the one standard deviation measurement errors and the summation is over 
the 42 A08 clusters), and so constrain parameters of given dark energy models. 

We construct a likelihood function L = e~ x2//2 , which depends on cosmological pa- 
rameters like Q m and those describing the dark energy model, as well as on the nuisance 
parameters. We marginalize over the nuisance parameters by multiplying the likelihood by 
the probability dis tribution function for the nuisance parameters and then integrating (e.g., 



Ganga et al.lll997l ). The resulting probability distribution function depends only on two vari- 
ables: Q m and a parameter p describing the dark energy model. In the ACDM case p is 
in XCDM it is uu x , and in 0CDM it is a. Since we consider only spatially-flat cosmologies for 
XCDM and 0CDM models, two parameters p and fl m completely describe the background 
evolution. 



2 These are discussed in depth by A08. Parameter A is the angular correction factor and is close to unity 
for all redshifts and cosmologies of interest (A08) so we take A = 1. 
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Since our initial likelihood function depends on 10 parameters in total (after setting A = 
1), to get to the two-dimensional probability distribution function we have to marginalize 
over 8 nuisance parameters (h, Q^, and the 6 parameters used to model the cluster gas 
mas s fraction). To re duce computational time we use the following statistics results (see, 
e.g.. lRiley et al.ll2002l . chap. 26). If two random variables a and b are independent and have 
probability distribution functions P a (x) and Pb(x), then variables c = ab, d = a/b, and 
/ = F(a) are also random with probability distribution functions 



Pc{x) 




P a (x')P b (x")5(x'x" - x)dx'dx" = I ^-P a (x')P b (x/x')dx', (2) 

' l x'\ 



Pd(x) = 




(3) 



PAx) 



dF-Hx) 



dx 



Pa(F-\x)), 



(4) 



where 5(x) is a Dirac delta function. For example, since g?a in Eq. (CQ) is inversely propor- 
tional to h, we define an auxiliary variable T = K^boiVt^h 2 ) / \fh and numerically compute a 
probability distribution function for it. This allows us to replace a five-dimensional integra- 
tion (over K, 7, bo, Qbh 2 , and h) by a one- dimensional integral over a new variable T and 
so reduce computational time significantly. 

Prior probability distribution functions for nuisance parameter are given in Table 1. 
The distribution functions for the 6 cluster gas mass fraction parameters are given A08, 
Table 4. Best fit values and confidence level constraints on dark energy parameters are 
sensitive to the assumed priors for the baryonic mass density and the Hubble parameter. 
Since different experiments give somewhat different estimates of these two, in our paper we 
use two prior sets to illustrate the differences. One set is from WMAP 3 year cumulative data 
JSpergel et alJl2007h . h = 0.73 ± 0.03 and Q h h 2 = 0.0223 ± 0. 0008, one standa r d deviation 
error s. The alternate set of priors we use is h = 0.68 ± 0.04 (IGott et al.ll2001t IChen et al. 



2003|) and VL h h 2 = 0.0205 ± 0.0018 ([Fields fc Sarkarll2006l ). one standard deviation errors. 



We define best fit values as pairs of values of cosmological parameters for which the 
likelihood function reaches its maximum. 1, 2, and 3 a confidence level contours are defined 
as sets of points in the parameter space with likelihood equal to e -2,30 / 2 , e -6,17 / 2 , and e~ n ' 8 / 2 
of the maximum value of the likelihood respectively. 
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3. Results and Discussion 



As a test we computed confidence level contours using our technique and the same 
prior set A08 used for ACDM and XCDM, Figs. 6 and 8 of A08. Our contours are in good 
agreement with those of A08. 

Figures 1 to 3 show cluster gas mass fraction confidence level contours and best fit 
values for the three dark energy models and the two sets of priors for fluh 2 and h given in 



Table 1. Compare d to the constraints d erived from the earlier lAUen et al.l (120041 ) cluster gas 



mass fraction data (jChen fc Ratrall2004l ). the difference between the contours corresponding 
to the two prior sets (for Vl^ti 2 and h) is much reduced. The new constraints are almost 
as restrictive as the ones derived from SNIa data and more constraining than those derived 
from angular size versus redshift data or Hubble parameter versus redshift data. 

Figure 1 shows constraints on ACDM. Q m is better constrained than Q\ and the results 
are in good qualitative accord with previous analyses. The best fit values are slightly away 
from a spatially-flat model. 

Figure 2 shows constraints on the XCDM parameterization. Again, the energy density 
of nonrelativistic matter is fairly well constrained while the equation of state parameter is 
less constrained. The best fit values are again not exactly on the u x = — 1 line, which 
corresponds to the spatially-flat ACDM case. 

Figure 3 shows constraints on the 0CDM model. Q m is better constrained than a. For 
both sets of priors there is an upper limit on a. The best fit values are on the a = line 
which corresponds to the spatially-flat ACDM case, but there is a big part of evolving dark 
energy parameter space that still is not ruled out. 



We acknowledge support from DOE grant DE-FG03-99EP41093 and INTAS grant 061000017- 
9258. 
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Table 1. Prior probability distribution functions for nuisance parameters 





Parameter 


Allowance 


Distribution 


WMAP prior VL h h 2 


Q h h 2 


0.0223 ±0.0008 


Gaussian 


WMAP prior h 


h 


0.73 ±0.03 


Gaussian 


Alternate prior Vl h h 2 


Vt h h 2 


0.0205 ±0.0018 


Gaussian 


Alternate prior h 


h 


0.68 ±0.04 


Gaussian 
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Fig. 1. — 1, 2, and 3a confidence level contours for the ACDM model. Solid lines (x for the 
best fit at fl m = 0.27, Q\ = 0.89) correspond to the WMAP prior and dotted lines (o for the 
best fit at Vl m = 0.27, f2 A = 0.88) correspond to the alternate prior. Dashed = 1 — 
line corresponds to the spatially-flat ACDM model. 
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Fig. 2. — 1, 2, and 3a confidence level contours for the XCDM parameterization. Solid lines 
(x for the best fit at Q m = 0.29, u x = —1.2) correspond to the WMAP prior and dotted 
lines (o for the best fit at f2 m = 0.27, uo x = —1.17) correspond to the alternate prior. Dashed 
u x = — 1 line corresponds to the spatially-flat ACDM model. 
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Q 

m 

Fig. 3. — 1, 2, and 3<x confidence level contours for the 0CDM model. Solid lines (x for the 
best fit at Q m = 0.27, a = 0.0) correspond to the WMAP prior and dotted lines (o for the 
best fit at f2 m = 0.26, a = 0.0) correspond to the alternate prior. The a = horizontal line 
corresponds to the spatially-flat ACDM model. 



